Plasma Formation Dynamics in Intense Laser-Droplet Interaction 
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Using three-dimensional, relativistic particle-in-cell simulations we study the ionization dynamics in intense 
laser-droplet interaction. Of particular interest is the laser intensity and frequency regime for which wavelength- 
sized targets are not fully ionized. We find that — despite the fact that the plasma created at the droplet surface is 
overdense — oscillating electric fields may penetrate into the droplet under a certain angle, ionize, and propagate 
in the just generated plasma. This focusing effect can be attributed to the local field enhancements at the 
droplet surface predicted by Mie theory. The penetration of the fields into the droplet leads to the formation of a 
highly inhomogeneous charge density distribution in the droplet interior, concentrated mostly in the polarization 
plane. We present a self- similar, exponential fit of the fractional ionization degree which depends only on a 
dimensionless combination of electric field amplitude, droplet radius, and plasma frequency, with only a weak 
dependence on the laser frequency. The strong effect of a tight laser focusing on the inhomogeneous charge 
distribution in the droplet is also discussed. 
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The typical scenario in intense laser-matter interaction is 
that within a fraction of a near-infrared laser cycle an over- 
dense plasma is generated (see, e.g., (T] [2)). Subsequently 
the laser is mainly reflected, and only within a thin skin layer 
higher charge states may be created. However, fast electrons 
may penetrate into the target and collisionally ionize the inte- 
rior further. Present-day short- wavelength free-electron laser 
(FEL) sources (at, e.g., DESY in Hamburg, LCLS in Stan- 
ford, or S ACLA in Japan) are used to probe via Thomson scat- 
tering (3H51 inhomogeneous, overdense plasmas generated in 
that way. The successful interpretation of the experimental 
results often relies on accompanying numerical simulations. 
Molecular dynamics is a powerful tool that is widely used to 
describe the ionization dynamics in small laser-driven clusters 
(SHE). However, for wavelength- sized targets such as droplets 
the influence of the target on the propagation of the incident 
electromagnetic wave needs to be taken into account self- 
consistently. This requires the solution of Maxwell's equa- 
tions together with the equations of motion for the charged 
particles. In the case of weakly coupled plasmas the problem 
can be reduced to the solution of the Vlasov-Maxwell system 
of equations, which is efficiently achieved using particle-in- 
cell (PIC) codes (2). 

In this Letter we are interested in the extreme nonlinear op- 
tics due to ionization of wavelength- sized droplets. As the 
laser field propagation is determined by the plasma density 
and the plasma is generated by ionization, the charge state 
and density distributions may be sensitive to the ionization 
dynamics. In fact, even the strongest present-day lasers can- 
not fully ionize heavier elements so that the assumption of a 
preformed plasma may be inadequate. Furthermore, the skin- 
effect may prevent the laser from penetrating into targets that 
turn overdense in the course of ionization so that, in general, 
a richly structured space and time-dependent charge distribu- 
tion develops (9). 

Spherical, wavelength- size, homogeneous dielectric or 
metal objects in plane- wave electromagnetic radiation clearly 



fall into the realm of standard Mie theory |[T0l and are of fun- 
damental importance in optics. In fact, Mie scattering theory 
was used to characterize rare gas clusters in recent FEL exper- 
iments on a shot-to-shot basis fTTl[T2ll . 

In our Letter, we study the "non-linear Mie domain" using 
a 3D relativistic PIC code with ionization included. We show 
that in a certain laser intensity regime the droplet target is nei- 
ther fully ionized nor are charges only created at the droplet 
surface. Instead, fields penetrate under a characteristic angle 
into the droplet, ionizing atoms in the polarization plane and 
triggering plasma waves which collide in a focal spot. We 
present results for the fractional ionization degree at various 
laser intensities, wavelengths, and densities, that turn out to 
follow an universal scaling law. Finally, we study the effect 
of laser focusing on the fractional ionization degree, clearly 
demonstrating that electromagnetic radiation may be "fed" 
into the droplet if the focus is off-axis, in that way matching 
the optimal angle of incidence. 

The ionization of an ion with charge state Z — 1 and ion- 
ization potential / due to the electric field E is implemented 
using the tunneling ionization rate formula fT3l 

(2E ch \ 2n ' k 2 h \E\ ( 2E ch \ 
with k = ^P^, E ch = n* = Z x [^l. Here, m is the 

h ' cn me ' y / ' 

electron mass and i# is the ionization potential of atomic hy- 
drogen. When an ionization event takes place a free electron 
at rest is created at the position of the ion. The energy needed 
for ionization is taken out of the field via an "ionization cur- 
rent" ji on parallel to the electric field at the ion location. The 
value of jion is such that ji on • E is the work spent on ioniza- 
tion per time step CUE]]. Energy conservation is taken care 
of during the whole process, i.e., if the remaining field energy 
after the ionization of particles in a cell is not sufficient for 
further ionization, this cell is not considered anymore during 
the current time step fT6l . 



2 




FIG. 1: (color online). Electron density in the beginning (a) and at the end (b) of the interaction with the laser pulse. He 2+ density (c) at the 
end of the interaction. All lengths are measured in units of the laser wavelength A, time is expressed in the unit of the laser pulse period Tl. 
Light colors correspond to low density, dark colors to high density. For better visualization of the droplet interior a quarter of it was cut-out. 
Laser and droplet parameters are given in the text. 



We start by presenting typical results from PIC simula- 
tions of the interaction of an intense, plane-wave laser pulse 
with an initially neutral He droplet. A spatial resolution 
of Ax = Ay = Az = A/100, 125 macro-ions and 250 
macro-electrons per cell were used. Absorbing boundary con- 
ditions for the fields and for the particles were employed 
in propagation direction, periodic ones for the other direc- 
tions. A linearly (in ^-direction) polarized 10-cycle sin 2 - 
laser pulse of carrier frequency ujo enters the numerical box 
through the boundary x = and propagates into the re- 
gion x > 0. The dimensionless vector potential amplitude 
a = \eA/mc\ = \eE/muoc\ was 0.5, corresponding to 
a laser intensity I ~ 5.2 x 10 17 W/cm 2 , the wavelength 
A = 2ttc/ujo was 800 nm (i.e. for a laser period Tl — 2.66 fs). 
The density of the 2R = 4A = 3.2/im diameter He-droplet 
was p — 0.14 g/cm -3 . If the droplet was completely preion- 
ized such a density would correspond to an electron density 
n e o = 24n cr , where n cr = 1.8 x 10 21 cm -3 is the critical 
density for 800-nm wavelength light. The droplet center was 
located at x = 4 A, y = z = 10 A. In the simulations presented 
in this Letter impact ionization was "switched-off". Test runs 
showed that for the intensities considered the effect of colli- 
sional ionization during the laser pulse is much smaller than 
that of field ionization. Moreover, for self-consistency colli- 
sional absorption should be taken into account as well if col- 
lisional ionization is. A recently introduced microscopic PIC 
code is capable of incorporating both, albeit so far only for 
smaller targets lfT7l[T8l . 

Figure[T]shows snapshots of the volume distribution of elec- 
tron and He 2+ densities. In the beginning the droplet is non- 
ionized and thus transparent for the leading part of the laser 
pulse. Later, as the field strength of the laser pulse increases 
in magnitude, ionization becomes more efficient, and an over- 
dense plasma is generated rapidly on the droplet surface as the 
pulse propagates over it, leading finally to almost full ioniza- 
tion of a thin surface layer. However, we also observe that a 
highly inhomogeneous density distribution inside the droplet 
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FIG. 2: (color online). Electron (bottom), He x + (middle) and He^ + 
(top) density in two perpendicular planes (k, E) (left) and (k, B) 
(right) at t = 7T L (a) and t = 12T L (b). All lengths are measured 
in units of laser wavelength A, densities in units of n e o. Laser and 
droplet parameters are given in the text. 



is formed, concentrated mostly in the polarization plane. In 
particular, there seems to be a focal spot (blue area in the po- 
larization plane in Fig. [lj> and[T]:). The fractional ionization 
de gree I r — 47r ^ 3 n — J n e {v) d r of the droplet at the end of 
the interaction is ~ 35%. 

He 1+ , He 2+ , and electron densities in the two perpendicu- 
lar planes (k, E) and (k, B) at times t = 7T L and t = 12T L 
are plotted in Fig. [2] Charge density builds up inside the 
droplet, starting from a certain region on the droplet surface, 
most clearly seen in the He 1+ -plot in Fig. [2^. At the later 
time in Fig. |2j) the charge density fronts merged already, cre- 
ating the focal spot of He 2+ density. Comparing the charge 
densities in the two planes shows that the ionization dynam- 
ics mainly takes place in the polarization plane (k, E). The 
corresponding distributions of the electric field components 
(longitudinal E x and in polarization direction of the incident 
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FIG. 3: (color online). Electric fields (E x top, E y bottom) in two 
perpendicular planes (k, E) (left) and (k, B) (right) at t = 7T L (a) 
and t = 12Tl (b). All lengths are measured in units of laser wave- 
length A, electric fields in units of the incident field amplitude. 

laser electric field E y ) are presented in Fig. [3] It is seen that 
an oscillating electric field penetrates into the droplet where 
in Fig. [2] the charge density is created. This is an interesting 
example for an electric field propagating in a plasma that is 
created by it in the first place. 

We attribute the fact that the field and ionization front dy- 
namics originate from a surface region under a certain angle 
6 > 7r/2 (with respect to k) to a local, time-dependent field 
enhancement on the droplet surface. In order to corroborate 
this statement, we show in Fig.[4]the radial electric field along 
the droplet surface in the polarization plane vs time and 6 as 
obtained from the PIC simulation (a) and according to Mie 
theory fTOl (b). Standard Mie theory is formulated for plane 
incident waves. However, as Mie theory is linear we syn- 
thesized our pulse via spectral decomposition and added the 
fields coherently. In the Mie simulation the droplet is assumed 
to be homogeneous and conducting, with a dielectric constant 
6 = 1 — n e o/n cr . Under such conditions Mie theory pre- 
dicts in the strongly overdense regime (where the skin depth 
is <5 e c/ujp <C R with uj p = ^e 2 n e o/m e eo the electron 
plasma frequency) that the electric field on the droplet surface 
is perpendicular to it. In Fig.|4]the time axis has been shifted 
such that t = corresponds to the moment when the max- 
imum of the incident laser pulse arrived at the droplet cen- 
ter. Both PIC and Mie result predict maxima of the electric 
field on the droplet surface for angles 9 /it g [0.4, 0.7]. The 
slight disagreement in the field distributions in forward direc- 
tion (small 0) is due to the fact that in the Mie calculation 
the droplet is assumed conducting (i.e., completely ionized) 
from the very beginning whereas in the PIC simulation there 
is not yet plasma at the rear side of the droplet (see Fig. [2]). 
The field enhancement predicted by Mie theory is in excel- 
lent agreement with the PIC results (~ 1.9 times the incident 
field). 

The propagation direction x of the field structures inside 
the droplet seen in Fig. is tilted with respect to k, lead- 
ing to the observed focusing effect. In order to interpret cor- 
rectly these structures, we project the field components in- 
side the upper half of the droplet onto %, E T = E x cos ip + 
E y sin ip, E n = —E x sin ip + E y cos (p, with ip the angle be- 
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FIG. 4: (color online). Radial electric field along the surface in the 
polarization plane vs angle and time, obtained from the PIC sim- 
ulation (a) and as predicted by Mie theory (b). The color indicates 
the electric field in units of the incident field amplitude E. The hori- 
zontal black lines indicate the angle at which the electric field at the 
droplet surface is highest. 



tween x an d k (see Fig|5p). The resulting field distributions 
for E T and E n are shown in Fig. [5]). As the values of E T 
are several times bigger than the values of E n , we identify 
the field structures as a longitudinal plasma wave. The nec- 
essary matching of the plasma wave to the electromagnetic 
field at the droplet surface results in the tilt of x with re- 
spect to k because the phase velocity of the plasma wave is 
smaller than c. The frequency spectra of the electric field at 
two points x = 3 A, y = 11.5 A and x = 4.5 A, y = 11 A 
inside the droplet are shown in the Fig. |5|:. They peak at the 
frequency approximately equal to the local plasma frequency, 
which may be estimated from the plot of the local density vs 
time in Fig.[5]l. 

Figure|6]collects all our simulation results for the final frac- 
tional ionization degree I r . Introducing the dimensionless pa- 
rameter r] = a/(R\/5e), it turns out that for all the various 
cluster sizes R > S e , densities uj 2 ^> uJq, laser intensities 
and wavelengths simulated, I r is well described by I r 1 — 
exp (—777). In our case of He we find 7 = 1560. Note that the 
species-dependence only enters via the ionization potentials 
/ in the tunneling ionization rate formula ([T]). Inserting the 

expression for the collisionless skin depth 5 e = cj — oj^ 

we obtain 77 = eE / (2TimR{uj 2 — cjq)) — eE / (27imRuj p v ), 
showing that there is only a weak dependence on the laser 
frequency. Indeed, for tunneling ionization the electric field 
amplitude matters, not the laser frequency. For small laser in- 
tensity, when only the thin skin layer on the droplet surface 
gets ionized, one expects I r = 4^3/3 ~ R~ x • This limiting 
case is contained in our formula. 

Before concluding we show results for a focused laser pulse 
with a focal spot size of 3 A. If the laser spot size is small 
enough or the droplet is situated out of the focus, I r will 
depend on the electric field values in regions on the surface 
where the coupling to the droplet interior is efficient. Re- 
sults for normal, oblique, and grazing incidence are shown 
in Fig. [7] I r varies from 1.8 to 5.7% upon moving the focal 
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FIG. 5: (color online), (a) The 7T< — 2.2 TV cv 

level of the electron 

density in the droplet central plane z = 10 for successive times 
{t = 5, 5.4, 6.1, 6.8, 7.4, 8.1, 8.8, 9.4, 10.1 Tl). The straight line 
indicates the projection direction %. (b) The distributions of the 
electric field component E T parallel and E n perpendicular to % at 
t = 8 Tl . (c) Frequency spectra of the electric field E T at the points 
n = (3, 11.5, 10) and r 2 = (4.5, 11, 10) inside the droplet, (d) 
I rtcr vs time at the points n and r 2 . All lengths are measured 
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FIG. 7: (color online). Electron (bottom), He ~ l ~ (middle), and He i ~ 
(top) densities at £ = 12. 5 Tl of a i? = 2A He-droplet for three 
different positions of the droplet central plane yd with respect to the 
laser focus at yo = 8: (a) yd = 8 (normal incidence), (b) yd = 7 
(oblique incidence), and (c) yd = 6 (grazing incidence). The inten- 
sity of the 10-cycle pulse was 5.2 x 10 17 W/cm 2 . The laser focal 
spot diameter was 3 A. All lengths are measured in units of the laser 
wavelength A, the densities in units of n e o- 
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FIG. 6: (color online). Universal curve I r = 1 — exp(— 777) (dotted) 
with 77 = a/(Rk/8l) = eE / (2nmR(u 2 p - l>q)) for the final frac- 
tional ionization degree of the He droplet after the interaction with 
a plane-wave laser pulse. The total laser pulse duration in all cases 
was 26 fs. The symbols for different a are indicated in the plot. For 
some of them, runs with different density or laser wavelength have 
been performed, as indicated directly by arrows. 



spot from the droplet central plane (normal incidence) to the 
droplet pole (grazing incidence). At normal incidence the am- 
plitude of the electric field in the range of optimal angles is 
too low to enhance the ionization process inside the droplet. 
As a result, an overdense plasma is generated at the front side 
of the droplet surface only. Under oblique and grazing inci- 
dence the field amplitude for the optimal range of angles is 
sufficient to trigger ionization inside the droplet. 



In summary, we showed that a strong near-infrared or 
optical laser pulse interacting with an initially neutral, 
wavelength- sized He-droplet may generate a highly inhomo- 
geneous charge density distribution. A plasma may not only 
be created within the skin layer at the surface. Instead, the 
electric field penetrates into the droplet interior for certain an- 
gles of incidence predicted by Mie theory. This causes ion- 
ization inside the droplet, mainly confined to the polarization 
plane. A particularly high abundance of He 2+ is observed 
where the ionization fronts and the trailing plasma waves col- 
lide. The fractional ionization degrees for various droplet and 
laser parameters are found to be in good agreement with a 
self- similar, exponential fit that is a function of a dimension- 
less combination of electric field amplitude, droplet radius and 
plasma frequency, with only a weak dependence on the laser 
frequency. We restricted ourselves to He but qualitatively sim- 
ilar ionization dynamics are expected for higher- Z material at 
higher laser intensities. The existence of certain angles of in- 
cidence for laser penetration also explains why we found a 
strong dependence of the fractional ionization on the droplet 
location in tightly focused laser pulses. 
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